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In this work, we present a fully quantum theory of the plasmonic nanolaser, based on the maser 
model. Theory can be applied both to the microlasers with high Q-factor cavities and plasmonic 
nanolasers. We show that the latter is essentially a thresholdless device. We obtain the statistics of 
quanta, spectrum of the plasmonic radiation and the second-order coherence degree in the steady 
state. The limits of the model applicability are discussed in detail. All results are compared with 
the full numerical simulation. 
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I. INTRODUCTION 

The modern story of the utilization of the nanoparticle- 
living plasmons started a decade ago, when the pos- 
sibility of the enhancement of surface plasmons in the 
nanoparticles (NPs) by optical gain in dielectric medium 
was predicted theoretically [T] and demonstrated in the 
experiments pQE]. 

It has been realized, that extremely tight field confine- 
ment in surface plasmon modes can be used to achieve the 
strong coupling regime with the emitters [3]. The metal 
NPs became a promising platform for implementation of 
the great variety of effects known in quantum optics and 
cavity QED [5HH] . In parallel, papers QH] claimed the 
lasing on the individual NPs, and a lot of work focused 
on the realization of smaller and faster sources of light 
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A metal NP adjacent to the active medium appeared as 
a basic unit cell of many applications of nanoplasmonics, 
such as metamaterials, nanosensing, optical logic, etc. 
Though, theoretical description of this object (which es- 
sentially constitutes a nanolaser, alternatively known as 
spaser) trails far behind the experimental progress. In 
ref.|17) we showed that semiclassical limit of quantum 
theory (see, e.g. [D> '20 ), which is eventually equivalent 
to the classical Maxwell-Bloch approach [3T] , is unadapt- 
able for the nanolaser. It corresponds to the so-called 
thermodynamic limit of the laser [22], which is not the 
case, since in the nano-sized resonant cavity most of the 
spontaneous radiation of the emitters goes directly to the 
laser mode. 

Several works, addressing operation of the small lasers, 
go beyond the classical approximation. Papers [531 121] 
are of great help for understanding the single-emitter 
laser operation. Unfortunately, these results cannot be 
directly transferred to the case of nanolaser, since the 
latter has relatively poor Q-factor and needs ss 10 2 — Id 3 



* ilfedorov@gmail.com 

t sarychev andrey@yahoo.com 



[9] [TB] emitters per NP to achieve sufficient gain. Treat- 
ing of the problem in Heizenberg picture is performed in 
[2"5] . This work perfectly illustrates where the features of 
quantum behavior should be neglected to obtain a sys- 
tem of equations, which is solvable analytically. In [17] . 
we have shown that the number of quanta populating the 
plasmon oscillator cannot exceed the ceiling, which is of 
order of tens, due to thermal limitations. That is, the 
nanolaser is truly a nonlinear noise device, where fluctu- 
ations play a crucial role. 

Here, for the first time, we make an attempt to de- 
scribe the plasmonic nanolaser in the formalism of the 
density matrix. We use a low-loss approximation, which 
holds for the plasmonic nanolasers with cavities start- 
ing from the FOM of > 10 2 , which is the case for many 
realizations [TTl - fT5I |2"o1 |2"7] . This allows one to employ 
the maser approximation and solve the master equation 
analytically. 

In the second section, we perform quantization of the 
plasmon a bit neater than it is done in rcf. lj. In chap- 
ter 3, we describe the model of the active medium and 
present the master equation. We solve it and derive the 
applicability conditions for our model. We explore the 
statistics of quanta and number of plasmons. Fourth 
chapter is devoted to the coherence properties of the plas- 
monic radiation. We compare our analytical model with 
the numeric simulation and discuss the results. 



II. PLASMON QUANTIZATION 

Since plasmons are bosons, they are expected to be 
quantized similar to electromagnetic waves. But energy 
of the plasmon is not just that of the associated electro- 
magnetic mode. Plasmons are electromechanical oscil- 
lations, and kinetic energy of electrons is an important 
contribution. 

In [T] , plasmon quantization is performed in a declara- 
tive manner, that is, without connection with the classi- 
cal laws of motion. Here, we do the same thing explicitly 
for the plasmon in the waveguide-like geometry, which 
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behavior is well known ([21], 02] )• We consider the sim- 
ple case of the horseshoe geometry, shown in the figjl] 
Charge and current distributions on its arms are found 
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FIG. 1. Schematic view of the horseshoe. White rows indicate 
current direction in the antisymmetric mode. Electric and 
magnetic fields are parallel to the X and Y axes. Depth of 
the horseshoe in the y direction is I. 

in the form 



J = J n cos(fc„z), 
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where J n and Q n here are some functions of time. Canon- 
ical variables of the system can be chosen as 
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where V — adl is cavity volume and oj p is plasma fre- 
quency of the metal. The Hamiltonian of n-th mode then 
reads 
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It is straightforward to check that the canonical equa- 
tions of motion with the Hamiltonian [3] reproduce the 
Ohms and the charge conservation laws. Now quantiza- 
tion can be performed in the usual way. We let q n and p n 
to be operators with the commutation rule [q n ,p n ] = i 
(h is taken to be 1 throughout the paper). Annihilation 
operator of the plasmon then is 
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{v n q n + ip n ) , [a n , a+] = 1, 



and Hamiltonian [3] takes the familiar form: 
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Operator of the amplitude of the electric field at the ar- 
bitrary point inside the cavity is then 
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where a; is a unit vector in the x direction. 



III. MODEL AND SOLUTION 

In this section, we describe and solve our model. Essen- 
tially, it is the well-known maser model (see, for example, 
[351 HH])- We shortly introduce it in the next subsection, 
and then solve it considering features of our problem. 

A. Master equation 

For interaction of the plasmonic mode with the active 
medium we employ the usual laser scheme. We treat 
every unit of the active medium (atom) as the four-level 
system (fig{2| which effectively represents both dyes and 
quantum dots (QDs). 
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FIG. 2. Level scheme of the active unit. Transition between 
levels a and b is resonant with the plasmonic mode. Pumping 
mechanism brings the system from the ground state \g) to \a) 
via the short-lived band c. 

In the pumping process, active atoms experience acts 
of excitation at random moments in time |30j . Gener- 
ally, such act brings the atom to some state in the upper 
energy band. After the fast and nonradiative diffusion 
down this band, the atom encounters a long-lived laser 
transition where it interacts with the plasmonic mode of 
the nanolaser. Then, mixture of states \a) and \b) decays 
to the ground level, and the atom leaves the interaction. 
The atom is then "worked out" , and is ready to reload. 

The lower state \b) is usually chosen to be fast decaying 
to the ground state. In that case, direct transition from 
I a) to \g) can be neglected, and population of the level b 
is zero. 

Lossless evolution of the plasmonic mode and the res- 
onant atom transition, in the rotating-wave approxima- 
tion, is governed by the Jaynes-Cummings Hamiltonian 



H = va + a + uiu z + g(a-a 
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Here, v and oj are frequencies of the horseshoe and atom 
working transition. Lowering operator of the atom tran- 
sition cr_ is connected with the inversion operator a z as 
follows: a z = [<T + ,cr_]. Coupling strength g depends on 
the position and orientation of the atom dipole moment 
II relative to the electric field [U 
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Dynamics of the full density matrix (DM) is described 
by the equation 



to 10 1 eV. For the good-confined modes value of g of 
order 10 -2 eV seems to be the upper limit. 



dp 
dt 
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where V — g(a-a + + <r+a) is Hamiltonian [7] in the in- 
teraction picture, and C is the Lindblad superoperator 
describing the dissipation in the atom and nanoresonator 
EH: 
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(10) 
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Here, the first term corresponds to the atomic transi- 
tions, denoted by thin straight numbered arrows in the 
fig. [2j The second term stands for the polarization decay 
of these transitions. The third term is responsible for the 
cavity losses. In the following, letters T, 7 and a without 
index apply to the working atom transition. 



Decay constants from 10 together with the coupling 
strength g determine the time scales of the dynamics. 
As noted earlier, loss rates of the second transition T 2 
and 72 are much larger than others and drop out. Then, 
only four scales r, 7, n and g are important. T is the 
rate of the spontaneous emission into the free space and 
ranges from 10 -8 eV for dyes up to the 10 -6 eV for 
QDs. Dephasing 7 is usually much faster and is about 
10~ 2 eV for QDs and 10" 1 eV for dyes. FOM of the 
plasmonic nanocavities ranges from 10 up to 10 3 and can 
be even higher at lower temperatures, or for the other 
types of cavities. That is, region of interest for k is up 



B. Solution 

In the present model, an excited atom interacts with 
the plasmonic mode individually, i.e. it does not know 
anything about the neighboring atoms. After this single 
interaction, atom can leave some of its initial energy in 
the nanolaser. If such events occur rarely, this enhance- 
ment is enough to sustain only, say, one quanta of energy 
in the resonator. More frequent "kicks" , providing energy 
to compensate larger losses, may lead to the laser action. 
It is a logical extension of the shots concept, which we 
used in our earlier work [T7| to describe the luminescence 
action in the nanolaser. On the language of the DM, the 
goal is to calculate the change of its elements as a result 
of a single excitation. 

In the present model, we consider the case when the 
plasmon lifetime k _1 is longer than time of interaction 
with the excited atom (we shall derive an explicit condi- 
tion later on). In that case, calculation of the "kick" can 
be performed analytically. One should solve equation [9] 
with an atom in initial pumped state and without the 
last term in|10| The plasmon dissipation enters only the 
coarse-grain master equation together with the constant 
flow of the kicks determined by the pumping rate. This 
approach is known as the maser approximation. 

Now let us turn on to the details of the calculation. For 
the plasmon oscillator, we use the basis of the occupation 
numbers. Elements of the full DM are indexed as p an ,pmj 
where a, (3 C a, b, g indicate the state of the atom. Mas- 
ter equation[9]then determines motion for these elements. 
In the maser approximation, it can be expressed in the 
form of a closed sets of equations: 
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Since T 2 is by far larger than other rates, equation for Pbn+i,bm+i can be excluded and one gets 
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dt 
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Equations 1 1 2 1 govern the internal dynamics of the kick. 
Initial value of the full DM is given by the outer prod- 
uct of the DM of the atom, prepared by the pumping 
mechanism, p®, and that of plasmonic field oscillator at 
the moment of atom excitation, 

P° P - P° = P° P ® Pa- In 
the simplest case, pumping brings the atom to the pure 
state \a). Initial conditions for the system 12 then are 

^ =0 

/i'n n rii *-* • 
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In the course of time, evolution [12] brings the atom to 
the ground state \g), so that p an ,am -> and p an ,bm+i -> 
0, Pbn+i,am 0. Final state of the plasmonic field can 
be obtained by tracing the ultimate DM over the atom 
variables: 
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Tr atom [pf mal ] =(g\pf" lal \g). (13) 
An expression for the kink of the DM then reads 
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n + m + 2s ' n + m + 2 + 2s 

where some indices are omitted for brevity: p n , m = 
(pT) , and we have introduced the constant s = 3s 

which is the saturation number of quanta [6]. Process 12 
has a lifetime 
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r + 7 - Rey (7 - T) z - Ag 2 (n + m + 2) 



which is analog of the atom transit time in the theory of 
micromasers. The whole interaction process is character- 
ized by the effective lifetime r e ff, which should be less 
than the plasmon relaxation time. This is the application 
condition of our model (f < 7): 
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where n is the average number of quanta in the oscillator. 

An important quantity which can be traced from here 
is the fracture of the spontaneous emission going to the 
laser mode. This value equals to the coefficient before 
Pn.m in the right-hand side of 14 if n = m = 0: 
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High value of (3 is distinctive feature of the microcavity 
lasers [221 In this respect, plasmonic nanolasers are 
beyond comparison in the optical range. As follows from 
the previous section, realistic numbers give s sa 10~ 4 and 
P as large as 0.9999. 

Now we can construct the large-scale master equation. 
In line with the maser approach, we sum up the atom 



kinks, expr 14 and add the plasmon dissipation term: 
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Here, R is the rate at which new excited atoms enter the 
interaction with the mode, that is, pumping rate. For the 
latter summand, we use the standard expression [28 : 



dt' 



p ) i oss = -n t h 2 (aa + p - 2a + pa + a + ap) 
- {nth + 1) g {a + ap - 2apa + + aa + p) , 



(19) 



where n t h is the number of quanta in thermal bath. 



Plasmon statistics 



First, let us determine the steady state of the main 



diagonal of the DM. From eq 18 we obtain a stationary 
condition for the diagonal elements — Pn^ dy state , 
which are probabilities of finding n quanta in the mode 
(for simplicity, we take nth = 0): 
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(20) 



Using the detailed balance condition, one finds the solu- 



tion: 



Pn — A 



T[n + s + 1] s^o 



(21) 



Here, T is gamma function and A is a normalization fac- 
tor. In case of strong coupling, when s — > 0, pff reduces 
to the Poisson distribution, regardless of the pumping 
rate. The total number of quanta in the mode is now 
straightforward: 
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FIG. 3. a) The average number of quanta [22] (blue line - theory; black dashed - numeric results), and condition 16 with values 
in electronvolts: k is yellow, t~}* is purple, b) statistics of quanta, at the pumpings marked by vertical dashed lines in the 
picture a (blue lines theory; black dots numeric calculation), s — 0.5 x 1CP 4 . Detailed paremeters: g = 10~ 2 eV, T = 10~ 7 eV, 
7 = W^eV, k = l(T 3 eV. 




10 100 1000 10 4 



a 




FIG. 4. The same as figjjjj but for s = 55.5 (g = 310 _5 eV, F = 10~ 5 eV, 7 = 10 _2 eV, « = 10 _7 eV). Horizontal dashed line in 
the fig. a is the threshold number of quanta. Again, condition |16| is fulfilled. 



where 1.F1 is the Kummer confluent hypergeometric func- 
tion. Thus, the strong coupling between the laser mode 
and the active media leads to the linear dependence be- 
tween n ss and R. Such behavior of laser is often referred 
as thresholdless, and is a commonplace in experiments 
on microlasers with high values of ft. 

Throughout the paper we shall illustrate our results 
using two sets of parameters. The first one stands for the 
real plasmonic nanolaser with Q-factor of 10 3 and strong 
coupling: g = 10~ 2 eV, which is the limit of s — > and 
P — > 1. The second one is closer to the nanolaser based 
on the photonic crystal or microsphere cavity. Our model 
covers these cases without any changes since the laser 
mode quantization for them leads to the same results of 
the section [n] Such cavities usually have much higher 
finesses and larger mode volumes, which results in the 
weaker coupling. We take j = 3x 10 _5 eV r and Q = 10 7 , 
which gives s = 55.5 and j3 — 0.017. 

Figure [3] illustrates the results [21] and [22] in the first 
case. In the fig[3^, the blue line shows the dependence 



between the average number of quanta and the pumping 
parameter a. It agrees perfectly with the result of the nu- 
meric simulation (black dashed line), which does not use 
the maser approximation. It is linear, as long as all the 
excitation energy goes to the mode, either in the spon- 
taneous or stimulated way. Still, care should be taken 
to satisfy the restriction [16] Yellow line indicates the 
present value of k and purple line stands for t~^. Since 
the purple line is above yellow one everywhere, condition 



16 is met and model is adequate for the parameters used. 

Blue lines in the picture [3]d show the statistics of plas- 
mons [2l] fo r three different pumpings, marked by dashed 
lines in |3^,. Black dots stand for the distributions ob- 
tained in the numeric calculation, which generally follow 
the theoretical prediction. We shall explore the differ- 
ences in more detail in discussion of the second-order 
coherence. 

Distributions are close to Poissonian in a good approx- 
imation, as expected in the limit ft —¥ 1: the first set of 
parameters gives ft « 0.99995. Note, that at the sec- 
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ond cross-section the number of quanta equals unity, and 
Pq s = Pf- These conitions express two different defini- 
tions of the laser threshold ( [251 1321 133] ) j which coincide 
in the limit in question. 

Figure [3] shows the same things for the second set of pa- 
rameters. Theoretical prediction of the number of quanta 
is surprisingly good, considering serious discrepancies in 
the statistics. Note, that the average number of quanta 
in the fig[4^i has a familiar "kink" at the threshold region, 
which indicates a relatively low value of ft, which amounts 
now to « 0.017. Limits in[2l1and[22lare no more relevant. 
As a result, distribution of quanta is different from the 
Poissonian, it is wider. Well above threshold, though, 
statistics is approaching that of the coherent state, as 
appears from the inset in fig]4]3. Section, marked by the 
second dashed line in [4^, is drawn in the point of max- 
imum slope of the blue line. At this point, number of 
quanta equals y/s (horizontal line in [4^,, as expected for 
classical good-cavity lasers [34] . 

It is easily seen, that all the DM elements, except those 
on the main diagonal, are zero in the steady state of |18| 
This state of the field is mixture of the pure states with 
the uniform distribution of phases. When elements of 
the main diagonal follow the Poissonian distribution, this 
mixture is often referred as the coherent state, though it 
is obviously not the case. Actually, DM of that kind just 
represents our state of knowledge of the field state with 
the total ignorance of its phase. Of course, it does not 
affect any of its coherence properties [35 . Interestingly, 
that the state of affairs changes, if atoms are pumped 
into the pure superposition of the working states. Then, 
phase of the field becomes partly predictable, which leads 
to the nonzero high-order diagonals in the DM [33] ■ 



IV. COHERENCE PROPERTIES OF THE 
PLASMONIC RADIATION 

What we are interested above all, are, of course, prop- 
erties of the outgoing radiation. In our model, it is rep- 
resented by the continuum of electromagnetic modes in 
the free space, constituting the reservoirs to which active 
atoms and the laser mode are coupled. Interaction with 
these modes is traced out and enters our model via the 
There, summands including the radiative de- 
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cay are those proportional to T\ and n. These constant 
actually account for the nonradiative decay also (ohmic 
losses in metal, for example) . The former, which is spon- 
taneous decay of the working atom transition into the 
free space, is almost entirely suppressed when /3 « 1. 
Therefore, we shall concentrate on the radiative part of 
the nanoresonator decay. 

It can be shown, that the corresponding DM of the 
electromagnetic oscillator in the far-field zone is exactly 
the same as for the radiating mode, up to the correspond- 
ing time delay and overall reduction of the intensity |28j . 
Accordingly, positive and negative frequency parts of the 
electric far-field are proportional to the plasmonic oper- 



ators a and a + . Methods, developed in the preceding 
sections, allow us to obtain the DM of interest; results 
on the quantum statistic thus hold for the outgoing ra- 
diation. That is, we have all instruments to explore the 
radiation of the plasmonic mode. 

The primary question to this section is: whether the 
radiation of the nanolaser is coherent? It is strongly con- 
nected with the definition of the laser threshold which 
became ambiguous as /3 is approaching unity. In this 
section, we shall address this issue from the framework 
of our theory. We concentrate first on the analysis of the 
first-order coherence, which will allow us to determine 
the spectrum and the linewidth of the plasmon radiation, 
and then turn to the study of the second-order effects. 

A. First-order coherence 

One of the widespread ideas about coherence is identi- 
fying it with the monochromaticity of radiation. This is 
a misleading view, although containing a grain of truth. 
In fact, monochromaticity does lead to the coherence of 
radiation, but only of the first order [35] and only for 
the stationary sources. Stationarity is obvious for the 
measurements in the CW regime but is questionable in 
the experiments dealing with the pico- and femto- second 
pulses. Here, we shall study the first-order coherence of 
the nanolaser radiation in the CW case, when the pump- 
ing rate can be thought as a constant. 

Both of these ideas rely on the concept of the first- 
order correlation function. Its normalized version, called 
the first-order degree of coherence, is defined as (28J [35] 

g l, )= (a + (t)a(t + T)) 
(a + a) 



(23) 



Here, operators a and a + mean the same as before, but 
now in Heisenberg representation are functions of time, 
taken at the same point in space. Angle brackets denote 
the quantum statistical average. Function [23] is valued 
between zero and unity. They correspond to the minimal 
and the maximal visibility of the interference fringes in 
the Young's experiment, which is the simplest procedure 
of observation the first-order coherence effect. Spectrum 
of the stationary field is the Fourier transform of this 
function l28l: 



S(u) =Re / drg^rje 1 - 



(24) 



That is, the challenge is to calculate the expression [24] 
with the two-time average in the numerator. As long as 
the motion of the DM, 
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be rewritten as follows 36 : 



is the Markov process, 23 



can 



Tr [all (t + r, t) (p(i)a 4 
(a+a) 



(25) 



where U is operator of the evolution 18 and p ss {t) is 



steady state of the DM. Multiplication by a + shifts the 
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values of zero diagonal to the first one. Function g 1 (r) zero steady state, as noted earlier. To obtain the decay 
is thus determined by evolution of the first diagonal of dynamics it is convenient to write its equation of motion 
the DM, which follows this same law [18] and leads to the in the following form: 



d_ 
dt 

where 



Pn,n+1 — R 



y/n(n + 1) (n + l)(n + 2) 



n + l/2 + s Pn ' n 1 n + 3/2 + s 



~Pn,n+l 



n = k [n + - - \Jn(n + 1) ] + R 



y/(n + l)(n + 2)p 



+ 3/2 - y/(n + l)(n + 2) 



n + 3/2 + s 



H"npn,n-\-l ) 

(26) 
(27) 



Eq. [26] without the last term has a steady state with 
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Pn,n+i r[n + 3/2 + s]' 



(28) 



The last term in [26] is a small disturbance relative to the 
rest of the system, and expression [27] is thus the decay 
rate of p n n +i, provided that statistics of these elements 
follows the quasi-stationary distribution [28] 

h , which has 



Initial point for the evolution 25 is p 



the nonzero elements only on the diagonal next to the 
main, whose distribution is 



Pn 



r[n + 3/2 + s]' 



(29) 



Form of this function is close to that of[28] and in a good 
approximation can be regarded as the same. That is, 
initial state for evolution U in [25] is already the steady 
one. Actually, it holds the same form [28] t hroughout the 
decay process since the disturbance in |26| is small. 

If n 3> 1, that is, in the limit of strong pumping, func- 
tion]^] is peaked around this large value, the only region 
where a substantial decay occurs; dependence of n in [27] 
then can be neglected and one obtains 



g 1 ^) = exp(-/vr) 

k R I L 



8n 8n 2 



(30) 
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As follows from [24] spectrum then has a Lorentzian shape 
with the traditional linewidth 30 Far more interesting 
is the spectrum in case of n < 1. In this domain, the 
approximation used in for the eq [30] goes wrong, as does 
the common model of the phase diffusion. Next, we shall 
dig a bit deeper and obtain an expression for g 1 which in 
valid for arbitrary values of n. 

In the evolution U (cq25), the decay with "statis- 
tics" [27] tends to distort the quasi-stationary distribution. 
This action is compensated at the price of small devia- 
tion of the actual distribution of elements p n , n +i from 
the [28] Therefore, the decay of the trace in [25] is indeed 
exponential for the p ss with any n. Its rate is then given 



by averaging of [27] over the steady state distribution [28] 
One obtains 



g 1 ^) = exp(-/i e //r), 



Veff 



1 X " 



(31) 



n=0 



That is, we expect the spectrum to have the Lorentzian 
shape with HWHM equal to p e ff- 

Sum in [3l] cannot be taken exactly in general. In the 
case a < 1, opposite to that of 30 analytical approxima- 
tion can also be obtained. Here, p n (27) can be fitted by 
the function starting from n = 1, and summation in 
[3T]then can be performed. Resulting expression is rather 
cumbersome, and here we show only its limit of small s: 



MS 



— > 
1 + a 



Ei(a) — In a + 7 e 

4~ 



(32) 



Here, Ei is exponential integral funciton and j e ~ 0.577 is 
the Euler's constant. Figure [5] shows these widths plotted 
versus a. Picture a stands for the case of small-s laser, 
with parameters as for fig[3] p e ff from 31 is plotted as a 
blue line. For comparison, we plotted the corresponding 
results obtained in the numeric solution of the master 
equation without the maser approximation (red dots). 

Essentially, in the fig[5[i one sees the linewidth pro- 
portional to the inverse number of quanta (recall that 
n ss oc a ), but limited from above by the cold-cavity 
linewidth, k/2. These two regimes, corresponding to the 
predominance of the spontaneous or stimulated energy 
transfer between active medium and the resonator, meet 
when n ss ~ 1. Expression |31| agrees with the numeric 
results perfectly, except the transition region, where the 
discrepancy reaches 20 percent. The purple line is the 
value of Pn i which approximates p e ff in the limit of 
strong pumping. The opposite case is covered by the 
approximation 



32 the green dashed line. 



Figure [5p is the same for the microlaser-like system, 
with parameters from the fig[4] We observe the kink- 
like element, resembling that in the plasmon number line 



a) b) 




FIG. 5. HWHM of the laser mode radiation, normalized to that of the free resonator. Blue lines: /i e //- Purple lines: Green 
lines: expr. |32| Red dots: numeric simulation. Plots a and b stand for the sets of parameters corresponding to plasmonic- and 
micro- lasers 



in figj4^. The purple line, or the "classical" linewidth, is 
perfect approximation after the kink, or classically, above 
the threshold. The slow-pumping approximation (green 
line) works even better than expected and holds up to 
a « s. Actually, we cover the linewidth analytically al- 
most in the whole range of a. 

To illustrate the logic of the spectra behavior clearer, 
we plot it versus a for several values of s on the same 
graph (figure [6]). 



the laser threshold, based on the line narrowing, fades 
away as we leave the thermodynamic limit. Line 3 
demonstrates a two orders of magnitude leap of the 
linewidth when the pumping rate doubles. No such 
thing is expected to happen in the plasmonic, thrcsh- 
oldless nanolaser. In that case, the only characteristic 
point which indicates the switch of the working regimes 
is a — n ss = 1. At that point, the linewidths roughly 
halves its initial value. The exact number is 0.53. 
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[3T| plotted for systems with 



FIG. 6. Linewidths 

10 _2 ,80, 10 3 (blue lines). Marked by numbers 1-3, respec- 
tively. Red dashed lines are proportional to the inverse num- 
bers of quanta. 

Red dashes, which essentially are the inverse numbers 
of quanta, lay precisely on the linewidth lines at large 
a. The departure starts when the linewidths begin to 
feel their ceiling, the cold-cavity value. This happens 
right after the kink-like intensity reduction (if any) oc- 
curs, when a ~ s (for large s). Origin of the lines in the 
fig |5] now became transparent: graph [5^, has the form of 
the first line in fig|6j while [5]d is closer to the second. 

Picture [6] also shows how the classical concept of 



B. Second-order coherence 

Now we turn to the analysis of the second-order coher- 
ence of the plasmon radiation. It is based on the features 
of the second-order correlation function 28 . 35 



9 2 (r) 



(a+{t)a+(t + r)a(t + r)a(t)) 
(a+a) 2 ' 



(33) 



Of greatest importance is its value for r = 0, which al- 
lows one to determine the nature of radiation considered. 
Signature of the macroscopic laser radiation well above 
threshold is g 2 (0) = 1, which characterizes the field in 
the m ixture of coherent states, discussed in the section 
IIIC As follows from eq 33 g 2 (0) is an intrinsic feature 



of the statistics of quanta, and can be readily calculated, 
given the distribution 21 Calculation results are shown 
in the figure [7J 

As appears from the graph a, our model misses essen- 
tial features of the statistics of quanta. Numeric result 
shows that the statistics is not Poissonian everywhere, as 
follows from the theory. Instead, it is super-Poissonian 
if a > 1 and sub-Poissonian if a < 1. Discrepancy 
manifests itself in the fig|3j), where the third dashed line 
stands nearly for the peak of the numeric line in[7^. Inter- 
estingly, that at the point a = 1, numeric line returns to 
the predicted value of unity. The reason why the theory 
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FIG. 7. Theoretical (blue lines) and numeric (black dashed) values of g 2 (0) versus a, for the first (a) and second (b) sets of 
parameters. 



fails to predict g 2 is that the condition 16 is not strong 
enough. As we have seen, it works well, if we are inter- 
ested in the spectrum of the radiation or the number of 
quanta. But inaccurate neglecting of the plasmon dissi- 
pation during the atom interaction time happens to be 
ruinous for the fine features of the statistics. To preserve 
them, sign < in |l"6| should really be replaced by <C. 

For the second set of parameters, condition [16] is satis- 
fied better, and consequently the figure [Tfj demonstrates 
a much better agreement. At weak pumpings, the theo- 
retical line starts from the value 



g 2 (0,a -> 0) = 2 



s + 1 



(34) 



For large s, it reduces to the value of two, which indicates 
the thermal state of light, in agreement with the results of 
the traditional laser theory. Kink in the g 2 (0) line occurs 
near the same value a as s, as it does in figures [4|5|6| At 
higher a, the value of g 2 (0) tends to unity, indicating the 
transition to the coherent state of light. 



V. CONCLUSION 

We have developed the density matrix model of the 
nanolaser, based on the maser approximation of the laser 
operation. We have shown that the plasmonic nanolaser 
is a truly thresholdless device, due to the extremely 
strong coupling of the plasmonic mode and the active 
medium, which result in values of j3, possibly as large 
as 0.9999. We obtain the analytical expressions for the 
number of quanta and the linewidth of the laser radia- 
tion, which are valid at any pumping rate and agree with 
the numeric simulations. Our theory also covers a broad 
range of the high-Q microcavity lasers. 

Plasmonic nanolasers, with respect to the Q-factor, are 
close to the edge of applicability of the maser model. 
The latter allows one to analyze the intensity of laser 
radiation and first-order coherence effects, but misses the 
fine features of the plasmon statistics, which determine 
degree of the second-order coherence. We conclude, that 
for the further theoretical DM-based investigation of the 
bad-cavity plasmonic nanolasers, one should go beyond 
the maser approximation. 
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